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O ' 

' We show that a Bose-Einstein condensate with a vortex lattice in a rotating anisotropic harmonic 

, potential exhibits a very low frequency scissors mode. The appearance of this mode is due to 

the SO(2) symmetry-breaking introduced by the vortex lattice and, as a consequence of this, its 
frequency tends to zero with the trap anisotropy e, with a generic dependence. We present 
analytical formulas giving the mode frequency in the low e limit and we show that the mode frequency 
for some class of vortex lattices can tend to zero as e or faster. We demonstrate that the standard 
classical hydrodynamics approach fails to reproduce this low frequency mode, because it does not 
' ' contain the discrete structure of the vortex lattice. 
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Since their first observations in 1995 the trapped gaseous atomic Bose-Einstein condensates have been the 

^ . subject of intense experimental Q and theoretical investigation. A fruitful line of studies has been the measurement 
' of the eigenmode frequencies of the Bose-Einstein condensate, which can be performed with a very high precision and 

a therefore constitutes a stringent test of the theory. Of particular interest in the context of superfluidity are the scissors 
I ' modes, first introduced in nuclear physics |5| and occurring in quantum gases stored in a weakly anisotropic harmonic 
I trap. These scissors modes have been studied theoretically 0] and observed experimentally 0] in a non-rotating gas, 
C5 ■ both in a condensate (where they have a large frequency even in the limit of a vanishing trap asymmetry) and in a 
O , classical gas (where their frequency tends to zero in a vanishing trap asymmetry limit). 

I Another fruitful line of studies has been the investigation of rotating Bose-Einstein condensates: above a critical 

rotation frequency of the anisotropic trap, it was observed experimentally that vortices enter the condensate and 
y—i , settle into a regular lattice 8, 9, 10] . In this paper, we study theoretically the scissors modes of a Bose-Einstein 
' condensate with several vortices present, a study motivated by the discovery of a low frequency scissors mode in 
, three-dimensional simulations in |lll |. In section we define the model and the proposed experiment, and we present 
^2 ' results from a numerical integration of the Gross-Pitaevskii equation: condensates with several vortices present have a 
' scissors mode with a very low frequency, unlike non-rotating condensates. We show in section^]that the equations of 
classical hydrodynamics, often used to describe condensates with a very large number of vortices 12], fail to reproduce 
our numerical experiment. Then, in the central part of our paper, section Ulll we interpret this mode in terms of a 
: rotational symmetry breaking due to the presence of the vortices, which gives rise for e = to a Goldstone mode of 
zero energy, becoming for finite e a low frequency scissors mode and this leads to an explanation for the failure of the 
classical hydrodynamics approximation. We then use Bogoliubov theory in section Hvl to obtain an upper bound on 
the frequency on the scissors mode, which reveals that two situations may occur: in what we call the non-degenerate 
case, the frequency of the scissors mode tends to zero as the square root of the trap anisotropy, and we show with 
perturbation theory that the upper bound gives the exact result in this limit. In the opposite degenerate case, the 
Q ■ frequency of the scissors mode tends to zero at least linearly in the trap anisotropy, and this we illustrate with a 
^ , numerical calculation of the mode frequency for a three- vortex configuration. We conclude in section IVl 



I. MODEL AND A NUMERICAL SIMULATION 

We consider a weakly interacting quasi two-dimensional (2D) Bose-Einstein condensate in an anisotropic harmonic 
trap rotating at the angular frequency 17 0|. In the rotating frame, the trapping potential is static and given by 

C/(x, y) = ^mcu^ [(1 - e)x^ + (1 + e)z/] (1) 
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FIG. 1: For an initial stationary condensate with vortices, numerical simulation of the excitation procedure of the scissors 
mode, by rotation of the trap by 10 degrees: density maps at successive times. Top row: 4- vortex configuration, with A'^ = 2500 
atoms and a trap rotation frequency Q — 0.7lo, resulting in ~ A.Ghio; the oscillation period is ~ 265ciJ~^ and the frames 
are separated by one quarter of a period. Bottom row: 10-vortex configuration, with A'' — 6000 atoms and a trap rotation 
frequency Q = 0.8w, resulting in ~ 5.9ftaj; the oscillation period is ~ 245ij;~^ and the frames are separated by one quarter 
of a period. The trap anisotropy is e = 0.025, the coupling constant is g — 0.0343/i'^/m. We have watched these oscillations 
without seeing significant damping over 7 periods. 



where m is the mass of an atom, e > is the trap anisotropy and uj is the mean oscillation frequency of the atoms. 
From now on, we shall always remain in the rotating frame. 

The condensate is initially in a stationary state of the Gross-Pitaevskii equation: 



2m 



A + 5|V'P + C/-f^iz 



(2) 



where g is the 2D coupling constant describing the atomic interactions |14|. /z is the chemical potential, ■0 is the 
condensate field normalized to the number of particles N and Lz = xpy — ypx is the angular momentum operator 
along z. In this paper, we shall concentrate on the case that the rotation frequency is large enough so that several 
vortices are present in the field tjj, forming a regular array. 

The standard procedure to excite the scissors mode is to rotate abruptly the trapping potential by a small angle 9 
and to keep it stationary afterwards. This is theoretically equivalent to abruptly rotating the field ■0 by the angle ^0 
while keeping the trap unperturbed: 



0,(^ = 0+) =e*''^^/V(i = 0-). 



(3) 



The subsequent evolution of is given by the time dependent Gross-Pitaevskii equation. 

We have solved this time dependent Gross-Pitaevskii equation numerically, with the FFT splitting techniques 
detailed in 0|, the initial state being obtained by the conjugate gradient method 0|. The results for a 4- vortex 
and a 10-vortex configurations are shown in figure ^ for e = 0.025, for an initial rotation of the trapping potential 
by an angle of 10 degrees: we see a small amplitude slow oscillation of the vortex lattice as a whole around the new 
axis of the trap, accompanied by weak internal oscillations of the lattice and of the condensate. This suggests that, 
in this low e limit, the excitation procedure mainly excites the scissors mode, and that this scissors mode has a very 
low frequency, much lower than the trap frequency oj. 



II. CLASSICAL HYDRODYNAMICS 



In the approximation of a coarse-grained vorticity, a condensate with many vortices is described by a density profile 
p{r) and a velocity field v(r) solving the classical hydrodynamics equations: 

(9f + V • grad) v = -— grad (U + pg) - 2^1 x v - U x (fl x r) (4) 

m 

dtp + divpv — 0. (5) 
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The first equation is Euler's equation in the rotating frame, including the Coriolis and centrifugal forces. The second 
one is the continuity equation. 

The velocity field v is the velocity field of the fluid in the rotating frame. In a stationary state, we set v = which 
amounts to assuming the solid body velocity field 17 x r in the lab frame. The corresponding stationary density profile 
is given by a quadratic ansatz: 

Po{x,y) = —{^J.- aix'^ - a2y'^). (6) 

g 

From Eq.Q one then finds 

a^,2^^[u;'ilTe)-n']. (7) 

What happens after the rotation of the density profile by a small angle? To answer this question analytically, we 
linearize the hydrodynamics equations around the steady state, setting p = pQ + Sp and v ~ Sv. 

dfSv = — — grad(5/9g — 2ri X (5v (8) 

m 

dtS p + div poS-v = 0. (9) 

At time t — — and, to first order in the rotation angle 9, Sp — 26mxy{ai — a2)/ g. The subsequent evolution 

is given by the time dependent polynomial ansatz: 

777 

5p = —[c{t)-Y-5A{t)Y\ (10) 
9 

5v = 5B{t)r (11) 

where 6A{t) is a time-dependent 2x2 symmetric matrix and 6B{t) is a time-dependent 2x2 general matrix. These 
matrices evolve according to 



SA = -ATrdB{t) - A5B{t) - 5B^{t)A (12) 

3 1^ 
-1 0, 



5B = 2 5A{t) + 2n\ \\ 5B{t) (13) 



where we have set A = ( 1 . The constant term evolves as c + pTiSB = 



02, 

First, one may look for eigenfrequencies of the system solved by SA and SB. For e = 0, this may be done analytically 
|l6j: one finds one mode of zero frequency, and 6 modes of non-zero frequencies: 



±2lu,± 



(2^2 _ ^2y/2 [(2^2 _ ^2)1/2 _ ^1 _ (^4) 



At this stage, the presence of a zero energy mode for vanishing anisotropy looks promising to explain the low frequency 
of the scissors mode. At weak but non-zero e, however, a numerical diagonalization of the resulting 7x7 matrix shows 
that the zero frequency is unchanged, whereas the others change very little. Analytically, one then easily finds the 
zero-frequency mode for arbitrary e: 

5A . 6B = (0 (15) 

Why is there such a zero-frequency mode in the classical hydrodynamics? We have found that this is because 
there a continuous branch of stationary solutions of the classical hydrodynamics equations parametrized by two real 
numbers a and (3: 

gp/m — p — a^x^ — a^y^ (17) 



where 



a+ = fli - f]/3-f ia/? (18) 
a_ = a2 -I- f^a + (19) 
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These real numbers are not independent since they have to satisfy 

al3{a + I3) + 2aai + 2/3a2 = 0. (20) 

This is a second degree equation for /3 at a given a so it can be solved explicitly, giving rise to two branches. One 
of them contains the v = stationary solution as a particular case, with a = f3 — 0; it terminates in a point where 
a+ = a- = 0. Each stationary solution on this branch has a zero-frequency mode. 

So what happens after the scissors mode excitation, in the classical hydrodynamics approximation? We have 
numerically integrated the linearized equations Eas. (|12ll3|l with the initial conditions specified above these equations 
(see figure ^ . We find a scissors mode oscillation, however at a large frequency close to the e = prediction 
v — [(2a;^ — fpy^^ — 51] , in disagreement with the simulations of the previous section [T^ . 
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FIG. 2: In the linearized classical hydrodynamics approximation, angle of the long axis of the condensate with the x axis of 
the trap, as a function of time, after the stationary density profile was turned abruptly by a small angle. Here e — 0.025 
and = 0.8tJ. No low frequency oscillation appears: the main oscillation frequency is ~ 0.37aj, in agreement with the lowest 
frequency of the e = limit given by Eg. 1114^ . 



III. ROTATIONAL SYMMETRY BREAKING: FROM GOLDSTONE MODE TO SCISSORS MODE 

Having failed to find the scissors mode of section^using the classical hydrodynamics approximation we now discard 
it and turn to a quantum treatment of the problem. Consider a stationary solution of the Gross-Pitaevskii equation 
with a vortex lattice present, in the case of an isotropic trap, e = 0. This solution clearly breaks the rotational 
symmetry S0(2) of the Hamiltonian. Since this is a continuous symmetry group, Goldstone's theorem guarantees us 
the existence of a degree of freedom behaving as a massive boson I^i particular, this implies the existence of 

a zero energy eigenmode of the condensate, corresponding to the rotation of the wave function in real space. The 
"mass" of the Goldstone boson can be shown to be 

M = dn J rL,^p. (2f) 

The variable 9 conjugate to the Goldstone momentum P has the physical meaning of being a rotation angle of the 
lattice with respect to a reference direction. Note that a quantum spreading of this angular variable takes place when 
time proceeds, as a consequence of the Goldstone Hamiltonian P^/2M. 

Now, for a non-zero value of the trap anisotropy, there is no SO (2) symmetry any longer in the Hamiltonian, so 
that the Goldstone mode is turned into a regular mode of the condensate, behaving as a harmonic oscillator with a 
finite frequency. This means that the variable 9 now experiences a potential. An intuitive estimate of this potential 
is obtained by taking the stationary solution, rotating it by an angle 9 and calculating its Gross-Pitaevskii energy in 
the presence of a trap anisotropy: 

n-^ + Vi9) (22) 
where we keep only the 9 dependent part of the energy: 



V{9) = -muj^e [{x^ - y^){l - cos 261) -I- 2{xy) sin 26*] 



(23) 
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where (...)= J ip* . . .i/j- We note in passing that the fact that ip is a local minimum of energy imposes that V{d) > 
for all rotation angle 9 so that 

(xy) = (24) 
{x^~y^) > 0. (25) 

Quadratising V{9) around 6* = leads to the prediction that the angle (9) oscillates when initially its value is different 
from zero. The resulting motion of the condensate is an oscillation of the lattice as a whole, i. e. a scissors mode with 
an angular frequency: 

scissors- Q^i^^^^ ■ y ) 

This results in a ^/e dependence of the mode frequency in the case when (a;^ — y^) does not vanish in the low e limit. 

It is now clear why classical hydrodynamics does not exhibit this low energy mode: since it does not break SO (2) 
symmetry, the corresponding Goldstonc mode is absent and so the low energy scissors mode does not appear. We 
note that, in cases where S0(2) symmetry breaking would occur in the hydrodynamics equations, the scissors mode 
frequency would be predicted to tend to zero, even if the motion of the particles is not quantized; this was indeed 
shown by an explicit solution of the superfluid hydrodynamic equation for a rotating vortex- free condensate exhibiting 
rotational symmetry breaking and the ^ dependence of the scissors mode frequency was seen experimentally | (T9l |. 

IV. ANALYTICAL RESULTS FROM BOGOLIUBOV THEORY 

A. Technicalities of Bogoliubov theory 

A systematic way of calculating the eigenfrequencies of the condensate is simply to linearize the time dependent 
Gross-Pitaevskii equation around a stationary solution, taking as unknown the vector {StpjSip*) where dtp is the 
deviation of the field from its stationary value. One then faces the following operator. 



(27) 



where the complex conjugate of an operator X is defined as (r|X*|r') — (r|X|r')*. The scissors mode we are looking 
for is an eigenstate of this operator, and the corresponding eigenvalue over h gives the frequency of the mode. 

At this stage, a technical problem appears, coming from the fact that the operator Cqp is not diagonalizable 
|2fl| |. This can be circumvented by using the number-conserving Bogoliubov theory, in which the eigenmodes of the 
condensate are the eigenvectors of the following operator: 

^ = ( n n* Kgp ( ^ ) (28) 



(29) 



Q*j '"^^ Q 

where Q projects orthogonally to the condensate wave function 



and Q* therefore projects orthogonally to 0*. We shall assume here that the operator C is diagonalizable for non-zero 
values of the trap anisotropy e. We shall see in the next subsection that it is in general not diagonalizable for e = 
when several vortices are present. 

One then introduces the eigenmodes of £, written as {uk,Vk), such that — l^fcP = 1; to each of these 

eigenmodes of eigenvalue corresponds an eigenmode of C of eigenvalue — e^, given by (w^, u^) [IJ. We shall assume 
here that the condensate wave function is a dynamically stable solution of the stationary Gross-Pitaevskii equation, 
so that all the are real. We shall also assume that the condensate wave function is a local minimum of the 
Gross-Pitaevskii energy functional (condition of thermodynamical metastability) so that all the Ck are positive. 
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B. An upper bound for the lowest energy Bogoliubov mode 



We can then use the following result. For a given deviation dip* of the condensate field from its stationary value, 
we form the vectors 



e. 
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-^Mg^^J.) (30) 



where q is an integer. As shown in the appendix^ we then have the following inequality: 



miuk efc < ^ (31) 



\eo,e_i, 



where the scalar product (,) is of signature (1,-1): 



u\ I u' 



vj\v'P= y (32) 

We apply this inequality, taking for SiJ; the deviation originating from the rotation of ip by an infinitesimal angle: 
expanding Eq.Q to first order in the rotation angle, we put 

eb = ( _,-n* r *!/,*) • (33) 



After lengthy calculations detailed in the appendix ^ we obtain the upper bound 

where 9 is the angle of polar coordinates in the x — y plane. 

In the limit of a small e, we assume that the scissors mode is the lowest energy Bogoliubov mode, as motivated in 
section Unl Ea. l|M|) then gives an upper bound to the scissors mode frequency. An explicit calculation of dgU shows 
that the upper bound Ea. H34|l coincides with the intuitive estimate Ea. H2()|l : 

dlU ^2muj'^e{(l)\x^ -y^\(l)). (35) 

This upper bound suggests two possibilities in the low e limit. In the first one, {4>\x'^ — y'^\4>) tends to a non-zero 
value, which implies that the scissors mode frequency vanishes at most as e^^'^ . In the second possibility, {4>\x'^ — y'^\4>) 
tends to zero; under the reasonable assumption of a {4>\x'^ — y'^\4>) vanishing linearly with e, this implies a scissors 
mode frequency vanishing at most as e. This second case we term 'degenerate'. 

We point out that the degenerate case contains all the cases where the stationary wave function ?/; at vanishing 
trap anisotropy has a discrete rotational symmetry with an angle 7 different from tt. To show this, we introduce the 
coordinates rotated by an angle 7, 

x' = a; cos 7 — ?/ sin 7 (36) 
y' — a;sin7 + y C0S7. (37) 

The symmetry of t/i implies that — y'^\(t)) — (0|a;'^ — y''^\(t>) and {(j)\xy\4)) — {4)\x' y' . By expansion of these 

identities, and assuming 1 — cos 27 ^ 0, one gets {4>\x'^ — y'^\(t)) — 0. 

An important question is to estimate the relative importance of the modes other than the scissors mode excited by 
the sudden rotation of the condensate. As shown in the appendix^ under the assumption that the scissors mode is 
the only mode of vanishing frequency when e — + 0, and anticipating some results of the next subsection, the weight of 
the initial excitation cq on the non-scissors modes behaves as 

'\? = 0{e^) (38) 



both in the degenerate and the non-degenerate cases. 
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C. Perturbation theory in e 

We now treat the anisotropic part of the trapping potential as a perturbation, dU — muJ^t{y^ — x'^)/2. The 
Bogohubov operator £, considered as a function of e, can be written as 

C^=Co + 5C (39) 

where 6C is a perturbation. Note that the exphcit expression of SC cannot be given easily, as it involves not only dU 
but also the effect of the first order change of the condensate wave function entering the mean field terms of C; but 
we shall not need such an explicit expression. 

Apart from eigenmodes of non-zero eigenfrequency, the operator Cq has a normal zero energy mode characterized 
by the vector e„ and an anomalous mode characterized by the vector Cq, in accordance with Goldstone's theorem. 
They are given by 



QoLztpo 
-QIL*M 

Qodfiipo 



(40) 
(41) 



where ipo is the stationary condensate field for e = and Qq is the projector on the e = condensate wavefunction 
(f>o- As shown in the appendix IXI one has indeed CoCn = and CoCa = en- Within the subspace generated by e„ and 
Ca, the operator Co has therefore the Jordan canonical form: 



-'Olsubspacc — \ 



(42) 



Note that the normal and the anomalous vectors are, up to a global factor, adjoint vectors for the modified scalar 
product Ea. (|32() . in the sense that (e„, Cq) = Mq where Mq is given in Eq.jSlJ for e = 0. 

As usual in first order perturbation theory, one takes the restriction of the perturbation SC to the subspace and 
one then diagonalizes it: 

r\ fSCnn 1 + SC 

Olsubspacc - ^g^^^ 5^^^ J . I4dj 

The notation SCan means than one takes the component of 6Cen onto the vector Cq. Since the adjoint vector of Ba is 
en/Mo one has 

dCan^ M^\er,,SCer,). (44) 
A similar notation is used for SCnn and SCaa- Forming the characteristic polynomial of the 2x2 matrix Ea. H43l) and 

1/2 

taking all the SC^s to be 0(e), one realizes that to leading order in e, the eigenvalues are ±SCan i that is they scale 

as if SCan + 0. 

Finally, we use the exact identity (e^. Ceo) — A^((/)|(9g[/)|0), proved in the appendix 1X1 noting that e„ and Co differ 
by terms of first order in e, expanding £ in £0 + SC and using the fact that (e, >Coen) — {SmCoe) — whatever the 
vector e, we get for the scissors mode angular frequency 

-Swisses = 5Can + Oie^) = + Oie^). (45) 

When (x^ — y^)o 7^ 0, this coincides with the upper bound Ea. (|34|l to leading order in e: the upper bound is then 
saturated for low e. 



D. Analytic expressions in the Thomas- Fermi limit for the non-degenerate case 

In the Thomas-Fermi limit (/i 3> fuv) there exist asymptotic functionals giving the energy of the vortex lattice as a 
function of the vortex positions in the isotropic jl4l | and anisotropic ,22] cases. Here we shall use these functionals to 
evaluate the frequency of the scissors mode in Eq. H45|) . 
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From the general expression (Eq.(2.12) in ^2^) we take the simpHfying assumption that the distances of the vortex 
cores to the trap center are much smahcr that the Thomas-Fermi radius of the condensate, to obtain 

AE = AEo + - e)x^ + (1 + e)y^^ i^-2n\ log^| + 2V2^^ 

-V2^nujriJ2^og(d,j/l) + 0{e'^) (46) 



where AE is the energy difference per atom between the configurations with and without vortices, rj — h/ ^2Nmg = 
hu / (ii\/2n) <C 1 (with /i being the Thomas- Fermi chemical potential of the non-rotating gas), {xi, yi) are the coordi- 
nates of the i*'' vortex core, dij is the distance between the vortex cores i and j and / is a length scale on the order of 
the Thomas-Fermi radius and independent of e p^ . AEq is a quantity independent of e and of the vortex positions: 
therefore we shall not need its explicit expression. It is assumed that ^ uj, and that 

/7r\ 1/2 

n>n„,^i^-j w|iog?7l (47) 

i.e. the rotation frequency is large enough to ensure that each vortex experiences an effective trapping potential close 
to the trap center. An important property of the simplified energy functional is that the positions of the vortices 
minimizing Ea. H46|l are universal quantities not depending on O when they are rescaled by the length A |25l |: 

2'm[Q — flm) 

Now we calculate (x^ — j/^) using the HcUman-Feynman theorem: 

±E^-^mcj'l{x'-y^M\' (49) 

where E is the energy per atom and (f> is the condensate wave function. We have also that 

^E{e = 0) = ^AE{e = 0) (50) 
de de 

since the vortex- free configuration is rotationally symmetric for e = 0. Wc thus obtain [2^ 

lim(0|(x2 - y'M = 2V2^r^^^-^ - vl (51) 

£—►0 LO ^ — ' 

i 

Using this formula we study vortex lattices with up to ten vortices. By considering first the case e = 0, we find 
that they are all degenerate except the case of two, nine and ten vortices, see figure for = 9 and figure Ob 
for riy = 10. Interestingly the e = 10-vortex configuration differs from the one of the full numerical simulation of 
figure^ By minimizing the simplified energy functional for a non-zero e, we find that the 10-vortex configuration of 
figurenceases indeed to be a local minimum of energy when e is smaller than 0.0182: this prevents us from applying 
Ea. (|45|l of perturbation theory to the calculation of the scissors mode frequency of the numerical simulations. The 
10-vortex configuration of figure|3|3 ceases to be a local minimum of energy when e is larger than 0.0173. The 10-vortex 
configuration minimizing the energy for 0.0173 < e < 0.0182 is a distorted configuration which breaks both the x and 
y reflection symmetries of the energy functional, but not the parity invariance, see figure Ot. 

The last point is to calculate the derivative of the mean angular momentum with respect to f2. In the Thomas- 
Fermi limit, for e = 0, and to first order in the squared distance of the vortex cores from the trap center, the angular 
momentum is given by l23 | : 

lim(L^) ~ huy — ^/2TTniuiriy^ xj + yf, (52) 

i 

where is the number of vortices. Using the rescaling by A to isolate the f2 dependence we get 

A (i^) (e ^ 0) = -V2^mu/2Ll^ ^ ^2 ^ ^2 (53) 

i 

where d}? /dO, = — A^/(ri— r2,„) < 0. In this way, I'scissors for low e can be calculated from Ea. H45() in the non-degenerate 
case analytically in terms of the equilibrium positions of the vortices: 



^scissors ~ 4e(r2 f^m) ^ ^2 _i_ „,2 " i^"^) 
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FIG. 3: Using the simplified energy functional Eg. 1461 we find the minimum energy configurations of nine vortices (a) and 
ten vortices (b) for e = 0. These configurations are non-degenerate: ~ Vi)/ IZii^i + Hi) — 0.03 for both of them. For 

0.0173 < e < 0.0182 a cross-over takes place between the 10-vortex configuration of (a) and the one of figure^ such a cross-over 
configuration is shown in (c), for e — 0.018. 



E. Numerical results in the degenerate case 



In the degenerate case we must go to higher order in perturbation theory to get the leading order prediction for 
the frequency of the scissors mode. A simpler alternative is to calculate numerically the frequency of this mode by 
iterating the operator Wfj^ starting with the initial guess eo defined in Ea. (|33|) 28] . The corresponding numerical 
results are presented in figure 0] in the non-degenerate case of two vortices, and in figure |S1 in the degenerate case 
of three vortices. In the non-degenerate case, t'scissors found to scale linearly with e for low e, as expected, and 
the corresponding slope is is very good agreement with the prediction Ea. H45|) : the agreement with the asymptotic 
formula Eq.|SlJ) is poor, which could be expected since the parameters are not deeply enough in the Thomas- Fermi 
regime |25|. In the degenerate case with three vortices, we find numerically that i^scissors scales as e'^/^ for low e, which 
is indeed compatible with Ea. (|34(l which leads to a scissor frequency upper bound scaling as e. The fact that a strictly 
higher exponent (3/2 > 1) is obtained shows that some cancellation happens in the next order of perturbation theory, 
may be due to the threefold symmetry of the vortex configuration. 
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FIG. 4: Frequency squared of the scissors mode as a function of the trap anisotropy e for a condensate with two vortices. 
N = 2500, VL = 0.6a; and g = 0.0343fi^/m resulting in a chemical potential /x ~ 4.9/itJ. The black circles are determined from 
the iteration of the operator . The red disks correspond to the numerical evaluation of the upper bound Ea. l|34|l . The solid 
line corresponds to perturbation theory: it is the numerical evaluation of the leading term in e of Eg. 1451 . which is linear in e 
with a non-zero slope since the two-vortex configuration is non-degenerate. The red dashed line is a guide to the eyes. 



CONCLUSION 



In this paper we have studied the problem of a non-classical scissors mode of a condensate containing a vortex 
lattice. In 2D simulations of the Gross-Pitaevskii equation we showed that, when such a condensate experiences a 
sudden rotation by a small angle of the anisotropic harmonic trapping potential, the orientation of its vortex lattice 
will undergo very low frequency oscillations of the scissors mode type. Motivated by these numerical results we 
searched for this mode using the well-known classical hydrodynamics approximation where the condensate density 
is a smooth function showing no sign of the presence of vortices and where their effect on the velocity field is taken 
into account through a coarse-grained vorticity. We showed that this approximation does not contain a low energy 
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FIG. 5: Frequency of the scissors mode as a function of the trap anisotropy e for a condensate with three vortices. A*' — 2500, 
Q = 0.7lo and g — 0.0343ft^/m resuhing in a chemical potential n ~ 4.6huj. The black circles are determined from the 
iteration of the operator . The red disks correspond to the numerical evaluation of the upper bound Eq.lj^U- The result 
of perturbation theory is not shown, since it vanishes for a degenerate configuration. The black line is given by the equation 

3/2 3/2 

V = 1.582e ' and results from a fit of the black disks by a power law e ' over the interval e < 0.01. The red dashed line is a 
guide to the eyes. 

scissors mode. We then were able to explain this discrepancy using the Gross-Pitaevskii equation that treats quantum 
mechanically the motion of the particles. In this case, the density profile of the vortex lattice breaks the rotational 
symmetry which naturally gives rise to a Goldstone mode in the limit of an isotropic trap and, for a finite anisotropy 
it becomes a low energy scissors mode. The existence of a scissors mode at a low value of the mode frequency is 
therefore a direct consequence of the discrete nature of the vortices in a condensate, which is itself a consequence of 
the quantization of the motion of the particles. 

We obtained quantitative predictions for the mode frequency using two separate methods. First we calculated an 
upper bound on the frequency of the mode using an inequality involving the Bogoliubov energy spectrum. Second, 
using perturbation theory in e, the anisotropy of the trapping harmonic potential, we showed that this inequality 
becomes an equality to leading order in e when the expectation value (a;^ — y'^)o taken in the unperturbed state does 
not vanish; in this case the frequency of the scissors mode tends to zero as e^^^, and we gave an analytic prediction 
for the coefficient in front of e^/^ in the Thomas-Fermi limit. However, in the cases where the expectation value does 
vanish (which we termed "degenerate"), the frequency will be at most linear in e. We have illustrated this using a 
three- vortex lattice where the frequency vanishes as e'^/^. Also, in the general case, we have shown that the relative 
weight of the non-scissors modes excited by a sudden infinitesimal rotation of the trap tends to zero as O(e^) so that 
the excitation procedure produces a pure scissors mode in the low trap anisotropy limit. 

Finally we point out that the existence of the non-classical scissors mode does not rely on the fact that we are dealing 
with a Bose-Einstein condensate per se. In particular, we expect, based on the very general arguments of section UTTl 
that this mode will equally be present in a Fermi superfluid containing a vortex lattice and that its frequency will 
tend to zero as ^/e in the non-degenerate case just like in its bosonic counterpart. 

Laboratoire Kastler Brosscl is a research unit of Ecole normale superieure and Universite Paris 6, associated to 
CNRS. This work is part of the research program on quantum gases of the Stichting voor Fundamenteel Onderzoek 
der Materie (FOM), which is financially supported by the Nederlandse Organisatie voor Wetenschappelijk Onderzoek 
(NWO). 

APPENDIX A: DERIVATION OF THE UPPER BOUND ON THE BOGOLIUBOV ENERGIES 

In this appendix, we prove several results used to control analytically the frequency and the excitation weight of 
the scissors mode. 

We first demonstrate the inequality Eq. H31|) . The key assumption is that the the condensate wavefunction is a local 
minimum of energy and does not break any continuous symmetry (which implies here a non-zero trap anisotropy 
e 7^ 0) . Then the Bogoliubov operator C of the number-conserving Bogoliubov theory is generically expected to be 
diagonalizable, with all eigenvalues being real and non-zero. We then expand eb on the eigenmodes of £: 
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where is the set of modes normahzed as J \ uk\'^ — \vk\'^ = 1, and which have here strictly positive energies Cfc > 
since the condensate wavefunction is a local minimum of energy |21(|. The b^s are complex numbers; the amplitudes 
on the modes of the JF_ family (of energies — Cfe) are simply b^. since eb is of the form (/. /*). Using the fact that for 
the modified scalar product Ea. l32|) . different eigenmodes are orthogonal, and each eigenmode has a 'norm squared' 
equal to +1 for and —1 for we get: 

(eo,ei> _ Efce^+ 2efc|6fc|2 ^^^^ 



(eo,e_i) E.,^,26-i|6fep- 

This is simply the expectation value of e^, with the positive weights ^|6fcp. Hence Ea. l31|l . 
Next, starting from cq given by Ea. H33|l . we have to calculate e_i and ei to obtain Ea. (|34l) . 

To calculate e_i, we take the derivative of the Gross-Pitaevskii equation Eq.Q with respect to the rotation 
frequency fl, which leads to 

Co. ( *■*) ^ (_*,) + . (A3) 

We have to get £ instead of £gp: since ip is normalized to the number of particles N, which does not depend on fl, 
we have 



iP*dn^ = Z7 (A4) 

where 7 is real, so that dntp = ijtp/N + Q{dQ,i^) where Q projects orthogonally to t/i. Using the fact that [ip^ —V'*) is 
in the kernel of Cqp, we obtain: 

^-1 = [ ^n*\f^i!*\\ (^^) 



which allows to get 

(eo,e_i) - / rL.iP. (A6) 



In passing, we note that the fact that all efc are positive implies that (eo, e_i) > so that the mean angular momentum 
is an increasing function of fi, a standard thermodynamic stability constraint for a system in contact with a reservoir 
of angular momentum 29] . 

We now proceed with the calculation of ei = Ce^. Since is hermitian, 

= Q{Lzij) - a^P (A7) 

where a is real. Since {4>, —"0*) is in the kernel of Cqp, we conclude that 



^Gpeo = i 





L! 



(A8) 



where [,] stands for the commutator. One calculates the commutator and then computes its action on (tp^—ip*): 
various simplifications occur so that 

f-h{deU)i,\ , . 

^''^^^ ^ [h{dgU)r ) ^ ' 

where 6 is the angle of polar coordinates in the x — y plane and where we used = —ihdg. ei results from this 
expression by projection. We finally obtain 



{eo,gi) = h' J ridm^p. (AlO) 

In what concerns the scissors mode frequency, the last point is to justify the identities CoCn — and CoCa = e„, 
where Cq is the e = limit of C, e„ is defined in Ea. H4U|l and Ca is defined in Ea. H41|l . First, one notes that, within a 
global factor i, e„ is the limit of eo for e — > 0. Then taking the limit e — > in Ea. ljA9|l . one immediately gets CoCn = 
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since dgU vanishes when the trapping potential U is rotationally symmetric. Second, one notes from Ea. ljASII that, 
within a global factor i, is the zero e limit of e_i. Then taking the zero e limit of the identity Ce^i — Cq leads to 

Finally we control the weight of the modes other than the scissors mode that are excited by the sudden infinitesimal 
rotation of the condensate: 



-♦non— scissors 1 1 2 ii \ ^ j / 

II =11 1^ h^^^ 



+ ^'U?* II' (All) 



where 1 1 . . . 1 1 is the usual £2 norm and fcs is the index of the scissors mode. The basic assumption is that the scissors 
mode is the only one with vanishing frequency in the e — s- limit. We rewrite Ea. (|A10|l as 

2efe|6fe|2 = 2mLu\ [ [x^ - y'^M'' . (A12) 



In the degenerate case the right-hand side is O(e^) so that each (positive) term efe|6fep is Oie^). For the normal modes, 
both e^^ and the mode functions u^-, Vk have a finite hmit for e — > 0, which proves | |g^ion^scissors| p _ Qi^^y 

In the non-degenerate case the same reasoning leads to a weight being 0(e). A better estimate can be obtained 
from 

(eo,ei)-e^Je'o,e_i) = ^ 2 (cfe - e^^e,) (A13) 
feej^+\{fe,} 

Then using the explicit expression for the scalar products in the left-hand side, see Ea. (|A6|) and Ea. ljA10|) . and using 
the result Ea. H45|) of the perturbative expansion for I'scissors = Cfc^/^, one realizes that the left-hand side of Ea. ljA13|l 
is O(e^), which leads to a weight on non-scissors modes being O(e^) as in the degenerate case. 
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